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1 Introduction 

The presence of hyperons in the inner core of a neutron star (NS) is still the subject of 
an open debate. The mechanism promoting the formation of heavier, distinguishable 
particles within dense interacting matter is quite intuitive at the non-relativistic level, 
and relies essentially on the Pauli principle. However, the energy and pressure reduc¬ 
tion due to the presence of hyperons in the medium seems to be incompatible with the 
most recent astronomical observations. In particular, most of non-relativistic models 
predict a softening of the equation of state (EOS) which is too large to sustain a 
neutron star of mass ~ 2M 0 as the ones that have been recently observed [1, 2], This 
apparent inconsistency between NS mass observations and theoretical calculations is 
a long standing problem known as hyperon puzzle. 

At present it is still too difficult to derive a first-principle effective Hamiltonian 
directly from LQCD results in matter with strangeness. Therefore it is still necessary 
to rely on the scarce experimental data to work out a realistic interaction between 
hyperons and nucleons in matter, and provide a reliable prediction for the EOS. We 
have been pursuing for a few years a program aiming to build up a phenomenological, 
realistic potential along the tracks of the Argonne/Illinois model. This work was 
started originally by Bodmer et al. [3], and continued in a series of papers by using 
variational Monte Carlo algorithms [4, 5, 6, 7, 8, 9, 10]. Our additional ingredient 
is the consistent use of auxiliary field diffusion Monte Carlo (AFDMC) [13, 14] 
calculations in oder to first assess the parameters of the interaction from measured 
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hyperon separation energies, and then to extrapolate the results to the case of neutron 
matter with strangeness. Other strange degrees of freedom might be relevant to the 
end of a correct determination of the EOS. However, in our approach we rely only 
on available experimental data. Essentially no measurements are available for E 
and 5 hypernuclei, and therefore we exclude at the moment these channels from our 
treatment. The same argument applies to two- and many-hyperon forces, which, at 
present, are substantially unknown from the experimental point of view. 

Most of our results have been recently published [15, 16, 17], and we refer to the 
original papers for all the detailed aspects of the algorithm and the actual calculations. 
In this proceeding we will briefly comment on a particular aspect of the hyperon- 
nucleon-nucleon force that nicely illustrates the difficulty in extracting the information 
on the Hamiltonian from experimental A separation energies. 


2 A phenomenological, realistic hyperon-nucleon in¬ 
teraction 


Within our non-relativistic many-body approach, A hypernuclei and A-neutron mat¬ 
ter are described in terms of pointlike nucleons and lambdas, with masses m jy and 
nr A, respectively, whose dynamics are dictated by the Hamiltonian: 


H nur — T + VnN ~ 
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Hhyp — H lmc + Ta + V AN + V\NN 


H nuc + ^ 2m A + ^ Vxi l ’ Xij ’ 

A A i A, i<j 


( 2 ) 


where A is the total number of baryons A = A/jv + A/"a, latin indices i,j = 1,..., ATn 
label nucleons, and the greek symbol A = 1,...,A/a is used for A particles. The 
nuclear potential includes two- and three-nucleon contributions while in the strange 
sector we adopt explicit AN and ANN interactions. 

In the non-strange sector we employ a simplified interaction in order to make the 
calculations feasible also for heavier hypernuclei. In particular we use the sum of 
the Argonne AV4’ interaction [11], plus the central repulsive term of the three-body 
Urbana IX potential [12], This choice provides a realistic description of energies, den¬ 
sities and radii. In Tab. 1 we report some examples of the binding energies obtained 
from AFDMC calculations for a set of closed shell nuclei. 

In the strange sector the AN interaction has been modeled with an Urbana-type 
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AV4’ 

AV4’TUIX C 

exp. 

diff. (%) 

4 He 

-32.83(5) 

-26.63(3) 

-28.296 

6 

16 0 

-180.1(4) 

-119.9(2) 

-127.619 

6 

40 Ca 

-597(3) 

-382.9(6) 

-342.051 

12 

48 Ca 

-645(3) 

-414.2(6) 

-416.001 

0.5 


Table 1: Energies of a few nuclei computed by AFDMC using the AV4’ nucleon- 
nucleon potential, and including a three-body repulsive contribution taken from the 
UIX potential. 


potential [18], consistent with the available A p scattering data 


V\i = v 0 (r Xi ) + ^v a T%(r Xi ) cr x - (Ti , 


(3) 


where Vq (r) = v c {r ) — vT%(r) is a central term. The terms v = (v s + 3vt)/4 and 
v„ = v s — v t are the spin-average and spin-dependent strengths, where v s and v t denote 
singlet- and triplet-state strengths, respectively. Note that both the spin-dependent 
and the central radial terms contain the usual regularized one pion exchange tensor 
operator T n (r) 
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3 
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where is the reduced pion mass 

m n o + 2 m n ± 

l-^7t ^ 


— ~ 1.4 fm . 

h"7T 


(5) 


All the parameters defining the AN potential can be found, for example, in Ref. [9]. 
The three-body potential v X ij can be conveniently decomposed in the 27r-exchange 
contributions +v^j S and a spin-dependent dispersive term v Xi j as follows: 


2tt,P ^P i x^ i 

V Xij = X iX , A X j f Tj • Tj , 


Aif = c sZ (r Xi ) Z (■ r Xj ) a, ■ r iX a t ■ f jX 


T i ' Tj , 


v \ij — W D T 2 (r X i) T 2 (r X j) 
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(7) 
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The function T n (r) is the same as in Eq. (4), while the X Xl and Z(r) are defined by 


where 


= Y w (r Xi ) cr x ■ a % + T n (r Xi ) S Xi , 
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AV?" 
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( 9 ) 


( 10 ) 


is the regularized Yukawa potential and S Xl is the usual tensor operator. The range 
of parameters Cp, Cs and Wo that have been used in our calculations can be found 
in Ref. [16]. 

AFDMC predictions for the A separation energy I3 X , defined as the difference in 
the binding energy of the core nucleus and the corresponding hypernucleus, are in 
good agreement with experimental data for hypernuclei up to 4 ®Ca and also for the 
A particle in different single particle states (see Fig. 1). 
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Figure 1: Measured [19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 
36, 37] and computed A separation energies as a function of /R 2 / 3 . Results for the 
A particle in different single particle states are also shown. 
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3 On the isospin dependence of the ANN interaction 


As shown in Ref. [17], the parametrization of the hyperon-nucleon potential yielding 
the best prediction for B\ provides a repulsion that is large enough to prevent the 
appearance of strange degrees of freedom in the range of densities found in the inner 
core of a NS. While this result might look as a possible solution of the hyperon puzzle, 
a closer look to the interaction suggests that the constraints coming from an analysis 
of the hypernuclear data might not be sufficient to provide an accurate extrapolation 
to stellar conditions. In the following we present a case study to effectively illustrates 
this point. 

The current version of the ANN potential does not depend on whether the two 
nucleons are in a singlet or a isospin triplet state. For symmetric hypernuclei the Pauli 
principle suppresses any strong contribution from the Ann or App channels. On the 
other hand, in neutron matter or in matter at /3-equilibrium the contribution of the 
isospin triplet channel might become quite relevant. In order to test the sensitivity 
of our AFDMC predictions for B\ on the strength of the isospin triplet component, 
we consider the sum w/T of the S and P wave 2-7T exchange terms: 


v x a T i 


t j = 


c 


^{x a ,X x ,\+C s Z( r M ) 


( 11 ) 


The operator r,; ■ Tj can be written in terms of the projectors on the triplet (T 
and singlet (T = 0) nucleon isospin channels: 


x — —3 P =u 4- P 

T j 6 r ij + 


T= 1 


The potential can be then rewritten as: 


Ti ■ Ti = -3 vlij P? 0 + C T v T Xij P? 1 


4 .. 


= 1 ) 
( 12 ) 

(13) 


where the additional parameter Ct gauges the strength and the sign of the isospin 
triplet contribution. For Ct = 1 the original parametrization of the three-body force 
is recovered. 

We have performed calculations for several hypernuclei ( 4 H, iHe, ;> He, ^O, 4 iCa, 
Ca) varying the Ct parameter in the range [—2; 3]. The results are reported in Fig. 2 
as the ratio of B\ (Ct) and B\ at Ct = 1, i.e. with the original parametrization. 

It can be observed that the sensitivity of the results on the value of Ct is not large 
over the whole interval —1.0 < Ct < 1.5, while strong deviations appear beyond this 
range. In general, outside of this range B\ tends to increase with respect to the orig¬ 
inal value. The only is exception is 4 ®Ca, where the sensitivity appears to be larger, 
and B\ tends to decrease with respect to the original value. Given the substantial 
asymmetry of this hypernucleus, we can infer that the isospin triplet channel could 
be properly constrained only by looking at the binding energy of strongly asymmetric 
hypernuclei. Extrapolating results to neutron matter without taking into account 
this feature of the interaction might in principle lead to misleading results. 
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Figure 2: A separation energies normalized with respect to the Ct = 1 case as 
a function of Ct- Grey bands represent the 2% and 5% variations of the ratio 
B\/B\(Ct = 1). Brown vertical arrows indicate the results for 49 Ca in the case 
of Ct = 2 and Ct = 3, outside the scale of the plot. 

4 Conclusions 

In the last years auxiliary field diffusion Monte Carlo has been used to assess the prop¬ 
erties of hypernuclear systems, from light- to medium-heavy hypernuclei and hyper¬ 
neutron matter. One of the main findings is the key role played by the three-body 
hyperon-nucleon-nucleon interaction in the determination of the hyperon separation 
energy of hypernuclei and as a possible solution to the hyperon puzzle. However, 
there are still aspects of the employed hypernuclear potential that remain to be care¬ 
fully investigated. For instance, we showed that the isospin dependence of the ANN 
force, which is crucial in determining the NS structure, is poorly constrained by the 
available experimental data. Future experiments on highly asymmetric hypernuclei 
such as 49 Ca, 9 J^Zr or even 209 Pb would pin down fundamental properties of the 
hyperon-nucleon forces. This would thereby allow for a substantial step forward in 
understanding the deep connections between the physics at the km scale typical of 
NSs and the properties of matter at the fm scale that can be efficiently explored in 
terrestrial experiments. 
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